Irrational mode locking in quasiperiodic systems 
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A model for ac-driven systems, based on the Tang-Wiesenfeld-Bak-Coppersmith-Littlewood au- 
tomaton for an elastic medium, exhibits mode-locked steps with frequencies that are irrational 
multiples of the drive frequency, when the pinning is spatially quasiperiodic. Detailed numerical 
evidence is presented for the large-system-size convergence of such a mode-locked step. The irra- 
tional mode locking is stable to small thermal noise and weak disorder. Continuous time models 
with irrational mode locking and possible experimental realizations are discussed. 
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Extended non-equilibrium systems subject to spatially 
uniform drives exhibit a vast range of spatiotemporal be- 
haviors, including uniform and static configurations, pat- 
terns periodic in space or time, and a great variety of 
spatially and temporally chaotic phenomena 1, 2\. In 
the cases where the local fields have a periodic attractor 
under a constant drive, the addition of an oscillating ex- 
ternal drive of period r can lead to mode locking: the 
local fields are translated by exactly p periods over the 
time qr, with fixed p,q € Z, for a finite range of ex- 
ternal drive amplitudes. This rational p/q mode-locked 
step is evident in the invariance of a dynamical variable 
(a velocity or current) to changes in the strength of the 
external driving force and is thus seen as a plateau in the 
average rate of change of the fields when plotted versus 
a control parameter. Rational mode locking has been 
seen in regular Josephson junction arrays driven by mi- 
crowave radiation (giant Shapiro steps) [3[ and a num- 
ber of translation-invariant fluid and granular systems 
4]. This mode locking can also be seen in systems with 
quenched disorder, e.g., experiments on [j| and models 
of [f| charge density waves (CDWs). A natural addition 
to these classes of spatially uniform or periodic systems 
and spatially random systems is the class of systems that 
have fixed quasiperiodic backgrounds. 

Spatially quasiperiodic systems have been extensively 
explored in a number of physical contexts. Quasiperi- 
odic solids are of course one example [7fl . The static and 
dynamic properties of mesoscopic quasiperiodic systems 
have also been a subject of study. In quasiperiodic wire 
arrays, the dependence of the normal-superconductor 
transition temperature on external field has cusps at irra- 
tional densities of the vortices [1, 0] . Similar complex be- 
havior has been seen in pinning structures which consist 
of quasiperiodic arrays of holes or mag netic (anti)dots 
added to a uniform superconductor d0|. Quasiperiodic 
arrays of colloidal particles can also now be constructed, 
using holographic optical traps [Til ] . 

We present here results on temporal mode locking 
in ac-driven two-dimensional systems with quasiperiodic 
structure. In addition to the typical rational mode lock- 
ing, our simulations provide clear evidence of mode- 



locked steps with p/q approaching an irrational value 
over a finite drive range in the large system limit. These 
steps are robust to small thermal noise and weak dis- 
tributed disorder. It would be of interest to search for 
these anomalous steps in experiments, where an output 
frequency would be irrationally related to the input fre- 
quency in the limit of large samples. 

For speed of computation, we carried out our initial 
simulations using a cellular automaton approximation of 
periodically driven overdamped systems. We first re- 
view its justification, based on arguments developed by 
Tang, Wiesenfeld, Bak, Coppersmith, and Littlewood 
(TWBCL) @,0,|l!,[l4|]. The automaton approximates 
a system with scalar degrees of freedom that represent 
particle displacements or coarse-grained coordinates for 
an elastic interface, such as the phase displacements in 
a charge density wave, flux vortex displacements along 
a drive direction, or fluid-fluid interfaces without over- 
hangs [lq . The essential features that are to be cap- 
tured are displacements that can be represented by an 
overdamped scalar field, local elastic interactions, spatial 
variation (such as impurity disorder that screens elastic 
interactions or a periodic background) that gives a natu- 
ral length scale, pinning by a potential that is periodic in 
the scalar field, and a periodic external force. These gen- 
eral physical features can be modeled by a continuous- 
time driven Frenkel-Kontorova (FK) model [l6j |. where 
N dynamical variables Xi (the coarse-grained interface 
or particle displacement indexed by lattice position i) 
evolve according to 

Xi = KA 2 x t + ti cos (kxi- fa) + F (t) , (1) 

where A 2 is the discrete Laplacian operator, K is the 
elastic constant, h! is the magnitude of the pinning 
potential, k determines the periodicity of the pinning, 
(3i G [0, 2tt) are pinning phases, and Fo(t) is a peri- 
odic driving force. Changing variables by substituting 
Ui = {kx t - A)/2tt, t -> 2nKt, h = h'/2nK, and 
F(t) = F a (t)/2nK gives 

ili = A 2 u. t + hcos(2irui) + F(t) + b { , (2) 

with bi = (27r) -1 A 2 fa in the range (— z, z), where z is the 
number of nearest neighbors [z = 4 for a square lattice). 
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Coppersmith [14| has shown that in the limit of large h 
and F(t) of the form F(t) = J2kez FS(t—kr) with ample 
relaxation time r between pulses, an automaton model 
is an appropriate approximation to these driven FK dy- 
namics. This claim is also supported by numerical work, 
which shows very similar mode-locking behavior for the 
continuous time model and the corresponding automaton 
@>[I3]- I n this limit, the dynamical variable Ui is near a 
minimum of the pinning potential, m ~ € Z, between 
drive pulses. One then defines an integer curvature vari- 
able Ci by 



A 2 n, + [bi + F — h\ 



(3) 



where [J is the floor function; this truncation of the net 
force approximates the selection of local minima rii by 
peaks in the pinning potential. The evolution of the Ci 
between subsequent pulses is the TWBCL automaton. 
At each time step, if the curvature exceeds the critical 
value z, the variable n, advances (due to the pulse and re- 
laxation) , redistributing the curvature ( "grains of sand" 
[Til ]) from a site i to neighboring sites j. Defining the 
toppling variable Ui(t) — 1 at all i where c,*(t) > z and 
Ui(t) = elsewhere, the update rule for a is 



Ci(*+1) = Ci(t)+A 2 Ui 



(4) 



In a given step of this automaton, the activity or toppling 
rate, corresponding to the velocity of the interface, is 
v t = N^ 1 Ui(t). Given periodic boundary conditions 
(BCs), the running dynamics of the system is fixe d by 
the initial conditions {c;(0)}, with = in Eq. (J3]) llal . 

In the Bak-Tang-Wiesenfeld sandpile version [13| 
(open BCs) and in disordered driven interfaces (reflect- 
ing BCs) , the initial c; are chosen from a Poissonian dis- 
tribution. We choose the initial Cj using Eq. ([3|) and 
quasiperiodic pinning 
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on a 2D square lattice of size N — L x x L y , i — (m, n), 
< m < L x , < n < L y , and frac(a;) — x — [x\ . The 
latter two terms in Eq. ([5]) maintain F — > — F symme- 
try. We choose the fractions P x /L x and P y /L y to be 
approximants to irrational numbers p x and p y , respec- 
tively. The approximants are found by truncating the 
continued fraction representation 
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We focus here on the case where p x = \[2 — 1 = 
(0, 2, 2,2,.. .) and p v = V3 - 1 = (0, 1, 2, 1,2,...). 

Given pinning {bi(0)} and drive F that fix a(t = 0), 
one characterization of the dynamics is the average ac- 
tivity rate v — (vt) (toppling rate or interface velocity). 
For any finite system there are only a finite number of 
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Figure 1: [color online] The activity staircase, velocity v vs. 
drive F, for the quasiperiodic automaton of size Lx x Ly — 
169 x 153 with pinning given by the rational approximants 
(70/169, 112/153) to (\/2 - 1,V3 - 1). Fractions with small 
denominators generally have the widest mode-locked steps, 
but the step at 86/1611 is an anomalous step which does not 
fit this pattern. Inset: Results for the FK-model, Eq. @, 
with h = 10 and the same size and quasiperiodic pinning; an 
anomalous step with similar v ~ 0.053 is seen. 



configurations, so there is at least one periodic attrac- 
tor (all attractors have the same v [l9(). We repeat 
the map, Eq. ([?]), until a periodic orbit is found and 
compute v(F) = M/(TN), where T is the period and 
M = ^2t—±^2iUi(t) is the number of topplings in one 
period. We find plateaus in v, i.e. mode-locked steps 
(Fig.[T]). Most of the wide steps seen are low- denominator 
rationals, but there is at least one extra "anomalous" step 
which has a large denominator for its width. 

We find such a step, similar in position, size, and ac- 
tivity, for all approximants to (p x , p y ) = (v2— 1, v3 — 1) 
with L x > 29 and L y > 56. The activity v(L x ,L y ) on 
this step is not fixed at a single rational (Table [TJ) ; the 
denominator and numerator for v(L x , L y ) grow with sys- 
tem size, but v appears to converge as (L x , L y ) — > (00, 00) 
(Fig. [2|) . This strongly implies a limiting anomalous step 
at an unusual activity rate. The width of the selected 
step also appears to converge (Fig. EJ). Scanning Table H] 
along the p x = \pl — 1 approximant direction (varying 
L x , fixed L y ), it can be seen that the numerators p a ^ for 
v(L x ,L y ) satisfy the recurrence relation 



Pa,b 



2Pa-l,b +Pa-2,b 



where a (b) is the number of terms used in the approx- 
imation of p x (respectively, p y ). The denominators q a ,b 
obey the same relation. Similarly, along the p y = V3 — 1 
direction (varying L y , fixed L x ) the recurrence relation 
that agrees with all data is 



Pa,b 



Pa, 6-1 + Pa, 6-2 if b is odd 
2pa,b-i + Pa,b-2 if b is even 



(8) 
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Table I: The system-size dependence of activity (velocity) 
v(L x ,L y ) at the anomalous step. The model is the TWBCL 
automaton with quasiperiodic pinning. The fractions P x /L x 
and Py/Ly, which determine the pinning, are approximants 
to \/2 — 1 and %/3 — 1, respectively. To highlight the recur- 
rence relation, Eqs. |7) IH). the fractions are not all written in 
lowest terms; in most cases, GCD(p, q) — GCD(L X , L y ). 



with the denominators q a .t also obeying the same rela- 
tion. Note that Eqs. (|7l8p are also the recurrence re- 
lations for the numerators and denominators in, respec- 
tively, the approximants P x /L x and P y /L y . These recur- 
rence relations suggest the temporal concatenation and 
synchronization of subsystems, but the precise construc- 
tion of such spatiotemporal behavior is unclear to us. 
Jointly solving Eqs. (|7|8|) using four velocities as seeds, 
we obtain the velocity v in the (L x , L y ) — > (oo, oo) limit 

7903 + 379^3 + 142\/2 - 181 -y/6 

v = , (9) 

156238 v ' 

which is an algebraic irrational number, as our conjecture 
for the mode-locking value in the limit of infinite system 
size for this choice of (p x ,p y ). 

We have found anomalous steps using other irrational 
pairs (p x ,p y ). For small system sizes, there are no steps 
that satisfy the recurrence relations, but as the system 
size increases, irrational candidates appear, with more 
than one candidate for a given (p x , p y ). It is even possi- 
ble that there are an infinite number of irrational steps 
with nonzero width in the infinite size limit for typical 
irrational pinning choices. It would be very surprising if 
other quasiperiodic pinning schemes (i.e., patterns for hi) 
or lattices did not also give irrational steps. We have not 
found a way to predict, for a given driving choice, the ac- 
tivity of any irrational mode-locked step or the number 
of such steps. One-dimensional automata with nearest 
neighbor rules exhibit few rational steps and in the sim- 
plest case should not exhibit irrational steps fl7j ]. 

Numerically more challenging versions of this model 
include simulations of the continuous-time FK model 
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Figure 2: [color online] The size dependence of (a) activity 
v(L x ,L y ), and (b) width, w(L x , L y ), of the anomalous step. 
Each curve is for fixed L y . The convergence of v and w as a 
function of system size strongly suggests that there is a well 
defined activity v and step width w in the limit (L x ,L y ) — > 
(oo, oo); we conjecture that v ^ Q, Eq. ©, and w 7^ 0. 



Eq. (J2). We have studied the FK model using quasiperi- 
odic pinning variables identical to those for the automa- 
ton. We find results similar to those for the automaton, 
with both rational and anomalous steps converging to 
fixed width in a sequence of system sizes (Fig. [TJ inset). 

For applications to experiment, it is natural to ask 
about the robustness of the irrational mode-locked step. 
To address this question we have studied the effects 
of thermal noise and quenched disorder. We model 
thermal noise by adding random topplings to the zero- 
temperature parallel dynamics of the system, Eq. ([4]). 
These random topplings occur at finite rate r] per time 
step at each site. This noise randomly increases or de- 
creases the displacement variable n, by one period. To 
implement this, we add backwards topplings, which take 
place when Cj < —z. The toppling variable Ui(t) used in 
Eq. (j4|) is modified: in addition to setting Ui(t) = 1 at 
all i where Cj(i) > z, one sets Ui(i) = — 1 at alii where 
Ci(t) < -z, with Ui(t) = for -z < c. t (t) < z. 

Automaton models wit h p eriod q > 2 are unstable to 
thermal noise in general [20(, but we find that for small 
thermal noise the change in the velocity is limited. Weak 
thermal noise, which leads to nucleation of out-of-phase 
regions (20| and subsequent motion of domain walls, de- 
stroys exact periodicity, but the anomalous steps can still 
be clearly identified (Fig. [3]). Numerical simulations and 
analysis of similar systems, such as depinning in CDWs 
with weak noise, also show small corrections to (v) on 
most of the step and rounding near small-denominator 
step edges [21 1. 

The effects of quenched disorder were included by mod- 
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Figure 3: [color online] The effect of thermal noise on the 
anomalous step for dimensionless noise values r\ = 10~ 5 , 
10" 4 , 10" 3 for system sizes defined in the key. The hori- 
zontal axis is the scaled force / which is linear in F with 
each zero-temperature step starting at / = and ending at 
/ = 1. For easier comparison, the zero-temperature activity 
v(L x ,L y ) has been subtracted. For each value of noise, there 
is similar rounding of the plateau for all system sizes. 



ifying the pinning variables bi, to model jitter in both 
pinning phases and coupling K. Weak distributed dis- 
order does not affect the activity v(L x ,L y ) significantly, 
with 97% of the step unchanged for 0.1% Gaussian noise 
in the bi with (L x ,L y ) — (169,153). However, irrational 
locking is sensitive to strong rare disorder: we find that 
randomizing the bi at even one site in large systems dras- 
tically reduces the width of the anomalous step. This 
is to be contrasted with Poissonian initial conditions, 
which are maximally random, and thus have steps that 
are not sensitive to strong disorder, though wide large- 
denominator steps are not present. 

Given the general applicability of the FK model and 
the generic features of the automaton approximation, our 
results suggest that irrational mode locking might be seen 
in a number of experimental systems. The most likely 
candidates are constructed mesoscopic systems, where 
thermal fluctuations are small enough and there is precise 
experimental control. As noted, quasiperiodic supercon- 
ducting wire arrays and Josephson-junction arrays ex- 
hibit interesting static features as a function of applied 
field; dynamic ac transport experiments might well be 
modeled by elastically or plastically interacting vortices 
|22l [23I ] that would realize the conditions necessary for 
irrational mode locking. Recently [ljj , quasiperiodic col- 
loidal particle arrays have been constructed using holo- 
graphic optical traps. Experiments for dc driven parti- 
cles in periodic arrays [241 ] suggest that experiments in 
a quasiperiodic background may be used to search for 
irrational mode locking in a collective system. 

We have studied a broadly applicable model of dy- 



namics in an extended system, where we have chosen 
a quasiperiodic spatial background. We find anomalous 
mode-locked steps that appear to converge to an irra- 
tional frequency relative to the drive frequency, in the 
limit of large system sizes. This step is stable to modest 
amounts of both thermal noise and quenched disorder. 
One avenue for comparison with experiment would be to 
study automaton models on quasiperiodic tilings (see 25] 
for the v — > limit). There are a number of experimen- 
tal systems that are closely related to this model and in 
which such irrational mode locking might be seen. 

We thank Ofer Biham for suggestions. This work was 
supported in part by NSF via DMR-0219292, 0606424. 
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